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We have studied properties of nonlinear waves in a matiiematical model of a predator-prey system with 

pursuit and evasion. We demonstrate a new type of propagating wave in this system. The mechanism of 

propagation of these waves essentially depends on the "taxis", represented by nonlinear "cross-diffusion" terms 

in the mathematical formulation. We have shown that the dependence of the velocity of wave propagation on 

the taxis has two distinct forms, "parabolic" and "linear". Transition from one form to the other correlates with 

changes in the shape of the wave profile. Dependence of the propagation velocity on diffusion in this system 

7-^ , differs from the square-root dependence typical of reaction-diffusion waves. We demonstrate also that, for 

^ . systems with negative and positive taxis, for example, pursuit and evasion, there typically exists a large region 

. ' in the parameter space, where the waves demonstrate quasisoliton interaction: colliding waves can penetrate 

G [ through each other, and waves can also reflect from impermeable boundaries. 

PACS numbers: 87.10.+e Biol and med physics: General theory and mathematical aspects 

Keywords: population taxis waves, population dynamics, wave-splitting 
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INTRODUCTION. 
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^~P • Population waves demonstrate an example of the reconciliation of random movement of individual organisms with ordered, 

deterministic behaviour of their communities. Study of the mechanisms of such reconciliation of disorder and order is important. 
To explain specific phenomena characteristic of various population communities, we have performed studies of pursuit-evasion 
waves in a mathematical model of a predator-prey system. Our results presented in ||l|] and here demonstrate that population 
taxis waves exhibit special properties, which provide insight not only into specific mechanisms of population dynamics, but also 
broaden the class of nonlinear waves to include quasi-solitons in reaction-diffusion-taxis systems. 
^ The spatio-temporal dynamics of a system of two interacting species are often described in terms of a reaction-diffusion 

k> i system of the type 

dZ d^Z 

-g^=9iP.Z) + D. — , (1) 

where P is the density of the prey population, Z is the density of the predator population, Di and D2 are their diffusion 
coefficients, the nonlinear functions f{P, Z) and g{P, Z) describe local dynamics, including growth and interaction of the 
species, whereas the diffusion terms describe their spread in space, e.g. resulting from individual random motions. However, one 
characteristic feature of living systems is their ability to react to changes of the environment, and to move towards, or away from, 
an environmental stimulus, behavior known as taxis. Examples are chemotaxis, phototaxis, thermotaxis and gyrotaxis Oyl]. 
Many models of spatial dynamics of populations take taxis into account, and its importance has been recognized in modelling 
various biological and ecological processes, including propagation of epidemics, bacterial population waves, aggregation in the 
cellular slime mold Dictyostelium discoideum, dynamics of planktonic communities and of insect populations ||2l|4|,|^. The 
existence of travelling waves, and also stationary spatially- inhomogeneous structures, in interacting populations with taxis has 
been demonstrated experimentally and theoretically idlTllsilQt llflllllll . 



Here we investigate a system of partial differential equations describing two spatially distributed populations in a "predator- 
prey" relationship with each other The spatial evolution is governed by three processes, positive taxis of predators up the gradient 
of prey (pursuit) and negative taxis of prey down the gradient of predators (evasion), yielding nonlinear "cross-diffusion" terms, 
and random motion of both species (diffusion). In this paper we consider the problem in one spatial dimension, x, using the 
equations 

dZ , , d^Z d dP 

^-.9(^,^) + ^^-'^+^^^- (2) 

The diffusion coefficients, D, for simplicity are considered constant, uniform and equal for both species, and ^ (-f §§) ^^'^ 
■^ (■^§j) are taxis terms |7|]; ft,_ is the coefficient of negative taxis of P on the gradient of Z, ft,+ is the coefficient of positive 
taxis of Z on the gradient of P. We choose as local kinetics functions f{P, Z) and g{P, Z) the Holling type III form used by 
Truscott and Brindley 11211 to describe the population dynamics of phytoplankton, P, and zooplankton, Z: 

f{P,Z) = f3P{l-P)-ZPy{p-' + iy^), 

g{P,Z) = ^ZP'-l{P'^ + y'-)-wZ. (3) 

It is known that these kinetics demonstrate "excitable" behavior, and the reaction-diffusion system Q has propagating soUtary 
wave solutions Il2l 11311 . We now show how inclusion of the taxis terms can alter the properties of such solutions. Though 
predator-prey systems, with one or both populations demonstrating "intelligent" taxis have been studied before, by means of 
individual-based Monte-Carlo simulations J14i llSll and by using partial differential equations lll6l llTl llsL Il9l |2(1 12111 . our 
objective here to isolate and identify the specific role of the taxis terms in creating novel behavior. We demonstrate a new type 
of propagating waves in the system (|5}. The mechanism of propagation of these waves essentially depends on the taxis of the 
species and is entirely different from waves in reaction-diffusion system. Unlike typical waves in reaction-diffusion systems, the 
pursuit-evasion taxis waves can penetrate through each other and reflect from impermeable boundaries. 

Behavior in the form of solitary propagating waves is typical for many spatially extended nonlinear dissipative systems. 
Solitary waves that remain unchanged after collision with each other are less typical and are known only for a rather narrow 
class of nonlinear dissipative media 12211 . In this respect, such waves are analogous to the solitons in conservative systems, 
whose study, as stable particle-like features of nonlinear systems, remain a key interdisciplinary topic of modern mathematical 
physics. In the present paper we demonstrate soliton-like behavior and also wave-splitting in a class of waves which can emerge 
in population dynamics models as a consequence of taxis. 

DETAILS OF THE MODEL AND NUMERICAL METHODS 

We have performed extensive numerical studies of equations (2,3). Three finite difference schemes were used, differing in 
their approximation of the taxis terms Lu = -S^u{x, t) — 4^^- In all cases, we approximated 

{Lu){x.i,tj) « [a{xi+i,tj) {S{xt+i,tj) - S{xi,tj)) - a{xi,tj) {S{xi,tj) - S{xi-i,tj))] /{Sx)^, 

where a{xi,tj) was given by one of the following; 
Scheme A: The central implicit scheme 12311 . 

a{xi, tj) = 0.5(w(a;i, tj+i) + u(xj_i, i^+i)) 

Scheme B: an "upwind" explicit scheme. For the positive taxis (pursuit), 

a[xi,tj) = u{xi-i,tj) if S{xi,tj) > S{xi-i,tj), 

a{xi,tj) = u{xi,tj) \f S{xi,tj) < S{xi-i,tj). 

For the negative taxis (evasion), 

a{xi,tj) = u{xi,tj) if S{xi,tj) > S{xi-i,tj), 

a{xi,tj) = u{xi-i,tj) if S{xi,tj) < S{xi-i,tj). 

Scheme C: an "upwind" implicit scheme. For the positive taxis, 

a{xi,tj) = u{xi-i,tj+i) if S{xi,tj+i) > S{xi-i,tj+i), 

a{xi,tj) = u{xi,tj+i) if S{xi,tj+i) < S{xi-i,tj+i) 

For the negative taxis, 

a{xi,tj) = u{xi,tj+i) if S{xi,tj+i) > S{x,^i,tj+i), 



a{xi,tj) = u{xi-i,tj+i) ifS{xi,tj+i) < S{xi-i,tj+i) 

So, the "upwind" schemes rather than use the mean between values of the variables subject to taxis at two neighbouring grid 
nodes as in the central scheme, select one or the other depending on the direction of taxis, i.e. sign of the gradient of the attractant 
(see e.g. 12411 for the discussion of upwind schemes). 

The majority of calculations were based on scheme C with discretization steps Sx = 0.1, St — 5 x 10^"^ for most figures 
or scheme B and C with Sx = 0.5 and St = 0.01 for large-scale parametric studies fig.|3land scheme fig. ^3 Selected control 
calculations used scheme B with smaller steps, down to Sx = 0.01, St = 4 x 10^^, and schemes A and C with Sx = 0.01, 
St = 10-3. 

Unless specified otherwise, we have calculated solutions to equations ( I2I3> with the following parameter values: ly = 0.07, 
(3 = 1, w = 0.004 for two different values of 7; 7 = 0.01, which allows propagation of purely diffusive waves, i.e. with 
h+ = h- = 0, D > 0, and 7 — 0.016, which prohibits propagation of such wave. 



CASE I: THE REACTION-DIFFUSION-TAXIS WAVES FOR 7 = 0.01 
The different mechanisms of wave propagation for 7 — 0.01, D — 0.04 
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FIG. 1: (a-f) The profiles of waves with D = 0.04 and different taxis coefficients h± (at h- = 0, h+ — 1 solitary wave solutions do not 
exist). The different shape of the profiles suggests that different propagation mechanisms are involved. Note the oscillatory onset of the pulse 
front when both h+ > 0, h^ > 0. (g,h) The oscillatory onset of the front of the pulse of (f), magnified. The horizontal lines are at the steady 
state levels, (Pq, Zq)- The theoretical value for the oscillation half-length is 1.256 . . .. 



Fig-fflshows the profiles of population waves in a purely reaction-diffusion case (a) and with addition of taxis terms (b-f). The 
taxis terms significantly change the shape of the profiles. The value of the pursuit coefficient /i+ has much more pronounced 
effect than the value of the evasion coefficient /i_. If only evasion (/i_ > 0) but no pursuit (/i+ = 0) is added, waves tend 
to retain the same shape as purely diffusive waves, with long and smooth plateaus. The addition of pursuit (/i+ > 0) adds 
distinctive features, in particular non-monotonic behavior of predators around the front and/or the back of the wave. 

Here we suggest a qualitative explanation of wave shape change in terms of the pursuit term (/i+) only. Ahead of the wave, 
the system is at its stable equilibrium. Consider the effect of a local increase of the prey density P above the equilibrium. The 
resulting flux of predators to the area, described by the taxis term with the coefficient /i+, will deplete the density of predators in 
surrounding areas, and the conditions of equilibrium will be violated. Decreased density of predators will temporarily encourage 
growth of prey, followed by influx of predators, and the same sequence of events occurs progressively at each point in the 
spatial (x) direction, constituting a travelling wave in the population pattern. Note that no diffusion of either prey or predators 
is required; the phenomenon requires only the presence of taxis terms in equation ^. The "excitable" character of the kinetics 
in (|3} leads to a strong magnification of the localized increase of the prey population, through the prey-escape mechanism (prey 
multiply faster than predators). This feature, of course, is also essential for solitary waves in purely diffusive systems. An 
unusual feature of taxis waves is the oscillatory character of the front, see fig. [flg-h). Since these oscillations have a small 
amplitude, they can be described by linearized theory. Indeed, in a steadily propagating wave with speed c, variables P and Z 
depend only on the combination ^ = .x — ct, and satisfy the "automodel" system 



d^P , d dZ dP 



0, 



d^Z , d dP dZ 



(4) 



The speed of the wave in fig-fflf) is c = 0.3535, and the steady-state values of the variables are Pq = 0.05703, Zq = 0.13480. 
With these parameters, a straightforward calculation gives solutions in the form (P, Z)(^) « (Pqj-^o) + Re [{Pi,Zi)e^^], 
\Pi,Zi\ ^ \Po, Zo\, with Ai,2 ~ 1.9925 + 2.5014i. This predicts the half-length of oscillations along <^ coordinate of 
7r/Im(Ai_2) ~ 1.256, in a good agreement with the observed shape, see fig. Qg,h). This means that these oscillations are 
not a numerical artifact. An important feature of solitary taxis waves is the uniformity of their shape, amplitude and speed: after 
a transient, these are the same, regardless of the details of the initial conditions. In this they are similar to reaction-diffusion 
excitation waves and different from solitons in conservative systems. 



Quaslsoliton interaction of pulses 






(a) h- =h+ =0 
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FIG. 2: Space-time density plots showing interaction of waves in (a) purely diffusive in (b) and (c) taxis cases. On all panels, the length of the 
interval L — 500, the time scale t G [0, 3000]. Black corresponds to P = 0.9, white to P = 0. 



We have found that the system ( I2l3t has a region of parameter space where solitary waves interact as if they were solitons. 
That is, they do not annihilate, as reaction-diffusion pulses usually do, but penetrate through, or reflect from, each other (since 
the waves are indistinguishable, these two terms describe the same result). Fig.|2shows results of simulations in an interval of 
finite length L with no-flux boundary conditions ^\x=o.l = and ^\x=o,l ~ 0. Two waves were initiated simultaneously, 
one at each end of the interval; the results are shown as density plots. In the purely diffusive case, panel (a), the waves annihilate 
at the collision. With the taxis terms included, panels (b) and (c), two types of quasisoliton interaction are possible: (b) - the 
waves penetrate through each other on collision, and are then reflected from the boundaries; (c) - after the waves penetrate 
through each other on collision or reflected from the boundaries, they split. 

Soliton-like interactions of sohtary waves have been observed in some reaction-diffusion systems with excitable kinetics, both 
in numerics I25[l26ll27ll28il29ll30ll and in experiments 13(113 III . Such interactions are always limited to narrow parameter ranges 
close to the boundaries between excitable and oscillatory (limit cycle) regimes of the reaction kinetics. In contrast, fig.|3]shows 
regions in the parameter space corresponding to different regimes of interaction and propagation of taxis waves described by 
equations ( I2I3> . Both the existence of steady propagating pulses and their ability to penetrate/reflect have a complex relationship 
with the kinetic and propagation parameters. However, it is quite clear that the ranges of parameters providing soliton-like 
behavior are not in any sense narrow, and do not require proximity to the oscillatory regime in the reaction kinetics. Although 
large enough /i_|- is typically sufficient for propagation of waves, quasisoliton behavior requires both /i+ and /i_ to be non- 
vanishing. 

Another qualitative difference of interaction of taxis waves from waves in reaction-diffusion excitable systems is that in the 
non-soliton regime (hollow circles on fig. |3} the waves demonstrate interference. After the impact (fig. 0t,b,c) they do not 
annihilate straight away, but penetrate through each other (fig.Sl,e) and only decay afterwards (fig.Ef). 
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FIG. 3: Parameter regions corresponding to different regimes of taxis waves (D = 0.04). (a)/3 = 1,7 = 0.01, w = 0.004, (b) h+ = h- — 1, 
w — 0.004, (c) (3 — 1, ^ — 0.01, ft- = 1, (d) /? = 1, 7 = 0.01, /i+ = 1. Solid circles: quasisolitons pulses with soliton interaction. Solid 
triangles: quasisolitons with the wave splitting. Hollow circles: stable propagation of pulses with nonsoliton interaction on collision. Hollow 
triangles: propagation of pulses with the wave spliting. Dots: no stable propagation of pulses. Crosses: oscillatory local kinetics. Shown are 
two-dimensional cross-sections in the parameter space, all through the point corresponding to the standard parameter values. 



Wave-splitting 



The phenomenon of spUtting of taxis waves can be observed both in the case of sohton-like interaction (fig.|3j solid triangles) 
and in the case of non-sohton interaction (fig.|3j hollow triangles). Fig. [S] shows space-time density plots of various splitting 
regimes, depending on h^, observed during propagation of one-dimensional taxis waves with non-permeable boundaries. 

With the increase of /i+ (fig. |5j)-c), the time interval between subsequent splitting events decreases. The backward waves 
splitting from the forward propagating wave may decay, or split themselves (fig.|5}. This chain of splitting events can lead to 
self-supporting, aperiodic or approximately periodic activity. A single act of splitting is shown in more detail on fig. |6l The 
splitting event is preceded by an oscillatory instability of the propagating wave, which can be seen as the growing oscillations of 
dynamic variables; see e.g. the evolution of the amplitude of the prey population, A(t) = max P(x, t), on fig.|6j). 

x£lO,L] 

So, on fig. Etc), we observe how a single solitary taxis wave leads to a self-supporting splitting-wave activity. To see better 
the properties of the advancing front of this activity, we have used the explicit decomposition of the system (|3 with respect to 
the translation group, as described in 13211 : 



dP 
'dt 
dZ_ 

'dt 

P{x.,t) 
dX 



— = fiP,Z)+D 



giP, z) 

c{t), 



D 



d^p 

d'^Z 



d_ dz__ ap 

dx dx dx ' 



d_ dP 

dx dx 



c{t) 



dZ_ 

dx ' 



(5) 
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FIG. 4: Dynamics of the nonsoliton interaction of taxis waves. D = 0.04, h- — 2,h+ — 0.85 
t t t 





(a)h- =2,h+=2.8 (b) ft_ = 2, /i+ = 2.9 (c) /i_ = 2, /i+ = 3 (d) /i_ = 2, /i+ = 3.1 

FIG. 5: Space-time density plot propagation and interactions of pulses with spliting. L = 600; t G [0, 3500] 



where P* has been chosen equal to 0.5, and a;* well within the simulation interval, to minimize the effect of boundaries. In other 
words, we consider Q in a moving frame of reference, whose origin has coordinate X{t) in the laboratory frame of reference. 
This frame of reference is chosen so that the position of the wave front defined as P{x,t) = P* is always at a fixed selected 
point a; = a;* in it. The advantage of this approach is that we can observe the evolution of the front for a long time interval, 
unrestricted by the medium size. The results of simulation of (|3 are shown on fig.0(this is a small part, actual simulation was 
much longer). We can see that each backfiring wave only slightly affects the propagation velocity of the first front. Also, in this 
way we have ensured that this advancing front continued for a very long time (up to < = 7000), thus presumably is stable, and 
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FIG. 6: Dynamics of wave splitting (h- = 2,h+ = 2.8). (a) Spacetime density plot (x- G [150,400], t G [200,650]). (b) A(t) 
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the process of back-splitting remained approximately periodic for all that time. 

Similarly looking splitting waves have been observed in excitable Belousov-Zhabotinskii chemical system ll33ll . where they 
were due to the medium anisotropy induced by atmospheric oxygen, in numerical simulations of FHN-type equations 13411 and 
in a three-component reaction-diffusion model describing blood clotting 13511 . 
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FIG. 7: The results of simulation of J5j (ft- = 2, h+ = 3). (a) The splitting waves in a moving frame (S,=x — X{t)) of reference attached to 
their advancing front, (b) The coordinate of the front as function of time. 



Propagation velocity of taxis waves 

Fig- m shows dependence of the propagation velocity of taxis waves on h+ at various fixed values of /i_. It also shows 
symbols of the corresponding interaction and propagation regimes. The graphs of the velocity have two branches, "parabolic" 
and "linear". Above a certain value of /i+, the linear branches are almost independent of /i_ (fig.|8^). The transition, with h^ 
increasing, from the parabolic to the linear branch correlates with a qualitative change of the shape of the wave. The dots on 
figllt graph mark selected values of h^^, the wave profiles for which are shown on fig.|9] The parabolic branch of this graph 
correspond to "double-hump" shape of the P{x) profile, and the linear branch corresponds to a "single-hump" shape. 

This result is confirmed by phase plane figures (fig.llOt. One can see that for ft._ = 2, the transition from the parabolic to the 
linear branch takes places at h+ = 1.2, which agrees with fig.|9jd). 



CASE II: THE REACTION-TAXIS WAVES FOR 7 = 0.016 



In the previous section, we studied properties of taxis waves when the local kinetics admitted purely diffusive waves in Q, 
i.e. for /i_ = /i+ = and D > Q. When the parameters of the local kinetics are changed, such propagation may become 
impossible. We have found that the taxis waves may still be possible in such cases. Specifically, we increased parameter 7, and 
found regimes for which waves can propagate only with the taxis terms, but not without them. 
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FIG. 8: Dependence of the propagation velocity of taxis waves on the pursuit coefficient h+ for various values of the evasion coefficient h- . 
Vertical bars and symbols on panels (b), (c) and (d) designate regions corresponding to the different propagation and interaction regimes. Solid 
circle: quasisolitons pulses with soliton interaction. Solid triangle: quasisolitons with the wave splitting. Hollow circles: stable propagation 
of pulses with nonsoliton interaction on collision. Hollow triangle: propagation of pulses with/by the wave splitting. Dot: there is no stable 
propagation of pulses. Dots on the graph of panel (c) correspond to the profiles shown on fig.|9| 
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FIG. 9: The profiles of taxis waves for h- = 2 and different h+, corresponding to various parts of the velocity graph, marked by dots on 
fig.!;). 



Note that the central difference scheme (A) only allows computations for D > 0, whereas our "upwind" schemes (B) and (C) 
can work for D = as well. 

Fig.lllh shows the dependence of the behaviour of taxis waves on /i_ and h^, for fixed 7 = 0.016 and D ~ 0.04. The point 
/i+ = /i_ = is in the region where stationary propagation of waves is not possible. Yet, for this local kinetics, in the (/i_ , h^) 
plane there are regions of both quasisoliton and non-soliton waves. Unlike the local kinetics of Part I, for this 7 = 0.016 we did 
not see any splitting waves in the explored parameter region. 

Now, we consider the case of no diffusion, D = 0, at 7 = 0.016 (fig. II lb). 

Dependence of the propagation velocity on h^ for this case is shown on fig. El Here we also see two distinct branches of the 
graph, and the transition between them correlates also with the transition from quasisoliton to non-soliton regimes of interaction. 

The correlation between parabolic-linear transition of the speed vs h^ dependence with double-single hump transition is 
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FIG. 10: Phase planes corresponding to various taxis waves, corresponding to the profiles shown on fig.|9| 

illustrated by fig.^]and corresponding phase planes on fig. 1141 

As in the previously considered case of 7 = 0.01 and D = 0.04, the non-soliton waves at 7 = 0.016 and D = demonstrate 
partial interference (fig.ll5>. After the collision (fig.llSh.b). the waves penetrate through each other (fig.H;-e), and then decay 
(fig.Ef-h, cf. fig. 13. 

MECHANISM OF QUASISOLITON INTERACTION 

We consider the details of the soUton-like interaction of taxis waves for the case D = 0. Fig. ^]shows a sequence of wave 
profiles during such an interaction. As noted above, a feature of taxis waves is the low level of predators ahead of the prey 
wave, as the predators are attracted backwards by the prey density gradient. This backward gradient of predators encourages 
the forward movement of prey (see Fig. I16h). The meeting of two prey waves creates a high peak of prey density (Fig. I16h-c). 
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(b) D = 
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FIG. 11: Parameter regions corresponding to different regimes of taxis waves for 7 = 0.016 





(h)h- 



(C)h- 



7 



FIG. 12: Dependence of propagation velocity of taxis waves on h+ for selected values of /i_, for 7 = 0.016 and D — 0. The vertical bars and 
symbols on the graphs designate different regimes of propagation and interaction. Solid circle: quasisolitons pulses with soliton interaction. 
Hollow circles: stable propagation of pulses with nonsoliton interaction on collision. Dot: there is no stable propagation of pulses. The dots 
on the graph of panel (b) correspond to wave profiles shown on fie. ll3l for various values of /i+. 

This higher local density of prey attracts predators, which abandon the margins of the collision zone (Fig. I16b.c). The local 
growth of predators causes escape of the prey from the center of the collision zone towards the margins abandoned by the 
predators. These events invert the gradients of the populations and re-create front structures on the margins of the collision zone 
(Fig. I16t.d). which then cause generation of two new, "reflected" taxis waves (Fig. I16b-h). which subsequently restore their 
normal amplitude. So, the phenomenon of reflection is driven by interaction of both the pursuit and evasion taxis terms, forming 
a positive feedback loop. The interplay and positive feedback between the two taxis terms can also be elucidated by considering 
a simple linear analogue of equations 0: 






which are obtained from (|5} by putting /1.2 
equivalent to a Schrodinger equation for u 






dZ _ d^P 

dt dx^ ' 



(6) 



/.fp 



I? = and removing nonlinearity from the taxis terms. System (|5) is 

1/2 
ih_ Z. This is consistent with both the oscillatory fronts of taxis waves 
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FIG. 13: . Profiles of taxis waves for h- — 7 and different h+, corresponding to the selected points shown by dots on the graph fig. ll2b V 



and their ability to reflect from each other The role of nonlinearities appears to be in selecting a unique amplitude and shape 
of propagating waves, and restricting, compared to (|6}, the relative values of h± that allow reflection. Adding diffusion in (|6} 
destroys propagating waves, but not necessarily in (|3 where its dissipative effect may be compensated by the nonlinear kinetics. 

As we have already noted, some properties of taxis waves are essentially different from those of solitary waves observed in 
excitable reaction-diffusion systems. This difference can be seen in Fig. E] which shows dependence of the wave propagation 
velocity on the diffusion coefficient D. This dependence is clearly different from oc D^^^ law obeyed by reaction-diffusion 
waves. There is a marked change of this dependence near the transition between annihilating and reflecting waves, which is yet 
another evidence of different mechanism of taxis waves, especially of quasisolitons. 

Conclusions. We have studied a spatially distributed predator-prey system of equations, in which, in addition to or instead 
of diffusion terms, we have included terms describing taxis of each species on the other's gradient; predators pursuing prey, 
and prey escaping predators. We have found that the taxis terms change the shape of the propagating waves and increase the 
propagation speed, which is an evidence of a different mechanism of propagation of these waves. In this change the major role 
is played by the pursuit terms. Also, the taxis terms can change the interaction between propagating waves, i.e. make them 
penetrate/reflect, rather than annihilate. For this effect, both pursuit and evasion terms are essential. We have also identified 
regimes of "wave-splitting", that can be observed both during normal propagation, and as a result of interaction of taxis waves. 

Quasisoliton interaction of taxis waves is not a purely mathematical exercise. In l36ll . soliton-like interaction of bacterial 
population taxis waves were observed in vitro, where colliding waves continued to propagate after collision without delay (see 
fig.Cll. 

However, apart from immediate application to predator-prey or similar systems, our results illustrate that inclusion of taxis 
terms can dramatically change the behaviour of the system. Reaction-diffusion systems are the most studied class of models 
describing waves and patterns in spatially distributed systems. Meanwhile mechanisms of spatial interaction in real system vary 
widely beyond simple diffusion; let us mention just the role of "cross-diffusion" in the propagation of forest boundary ISTIl . in 
chemical and biological pattern formation L38..39. ■4(1 . turbulence-shear flow interaction in plasmas 12 111 , and cross-diffusion 
interaction between displacement and velocity in the Burridge-Knopoff model of tectonic slips 14111 . We hope that our present 
results may stimulate more attention into properties of dissipative nonlinear PDF systems with spatial terms beyond simple 
diffusion term, as these constitute a new class of models with significantly different properties. 

This study was supported in part by FPSRC grant GR/S08664/01 (UK) and by RFBR grant 03-01-00673 (Russia). 
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FIG. 15: Dynamics of the nonsoliton interaction for 7 = 0.016 and D = 0, /i_ = 5, /i+ = 2, time interval =5. 
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